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V) Abstract 

The coma of Comet 103P/Hartley 2 has a significant population of large particles observed 

p. as point sources in images taken by the Deep Impact spacecraft. We measure their spatial 

and flux distributions, and attempt to constrain their composition. The flux distribution of 

these particles implies a very steep size distribution with power-law slopes ranging from —6.6 

^h to —4.7. The radii of the particles extend up to 20 cm, and perhaps up to 2 m, but their 

exact sizes depend on their unknown light scattering properties. We consider two cases: 

bright icy material, and dark dusty material. The icy case better describes the particles 

if water sublimation from the particles causes a significant rocket force, which we propose 

as the best method to account for the observed spatial distribution. Solar radiation is a 

plausible alternative, but only if the particles are very low density aggregates. If we treat 

■^- the particles as mini-nuclei, we estimate they account for < 16 — 80% of the comet's total 

water production rate (within 20.6 km). Dark dusty particles, however, are not favored based 

on mass arguments. The water production rate from bright icy particles is constrained with 

_il an upper limit of 0.1 to 0.5% of the total water production rate of the comet. If indeed 

icy with a high albedo, these particles do not appear to account for the comet's large water 

production rate. 
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S_i 1. Introduction 



Comet 103P/Hartley 2 (hereafter, 103P or Hartley 2) is a hyperactive comet. The ac- 
tivities of most comets are consistent with surfaces that are over 90% inert, i.e., activity is 
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restricted to a few localized sources. In contrast, the gas production rates of the hyperactive 
comets suggest activity over ~ 100% of their surfaces or more. Comet Hartley 2's peak 
water production rate during the 1997 perihelion passage was Qh 2 o ~ 3 x 10 28 molec s _1 
(Combi et al. |2011). The sublimation rate of an isothermal nucleus with a 4% Bond albedo 



at 1.03 AU (the 1997 perihelion distance of Hartley 2) is 2.9 x 10 17 molec cm 2 s 1 (|Cowan 



and A'Hearn|l979 ). Thus, approximately 10 km 2 of surface area would be required to match 
the measured water production in 1997. However, constraints on the size of the nucleus by 



Groussin et al. (2004) and Lisse et al. (2009) suggested the total surface area was near 4- 



6 km 2 , yielding an active fraction (area based on Q// 2 o/ area of the nucleus) near 1.7-2.5. In 



order to account for the high active fraction, Lisse et al. (2009) proposed that the coma con 



tained a population of icy grains, which increased the surface area available for sublimation, 
and therefore allowed for a relatively high water production rate. 

The Deep Impact spacecraft flew by Comet 103P on 4 November 2010. The minimum 
fryby distance was 694 km, occurring at 13:59:47.31 UTC with a relative speed of 12.3 km s _1 



(A'Hearn et al. 2011 ). Rather than just being the fifth comet nucleus to be imaged up close 



by a spacecraft, Comet Hartley 2 might also be considered an archetype of the hyperactive 
comets. The flyby images verified the small surface area of the nucleus (5.24 km 2 ; Thomas 



et al.|this issue ), and Deep Impact's IR spectrometer revealed a coma of icy grains (A'Hearn 
et al. 2011). Knight and Schleicher (this issue) discovered tail of OH and NH, evident on 



10 km scales, which they conclude are derived from small grains of ice that were accelerated 



down the tail before completely sublimating. Bonev et al. (this issue) observed enhanced 



heating of the water gas 20-150 km down the tail, which they also attribute to sublimating 



ice grains in the tail. The icy coma hypothesis of Lisse et al. (2009) seems to be qualitatively 
valid. What remains is a quantitative verification that the ice grains in the coma are a 
significant source of water gas. 



A'Hearn et al. (2011) also discovered bright point sources distributed around the nucleus 



of Hartley 2 in visible wavelength images. The fluxes of these sources are consistent with a 



population of objects centimeter-sized or larger. A'Hearn et al. (2011 ) could not determine if 



these large particles were icy or refractory. In the following paper, we will attempt to discern 
the nature of these point sources and if they might be the cause of Hartley 2's hyperactivity. 

The point sources around Hartley 2 are found in visible wavelength images taken with 
both the High Resolution Instrument (HRI) and the Medium Resolution Instrument (MRI). 
Their spatial distributions and densities vary, but they are not background stars; celestial 
sources are significantly streaked in images taken while the spacecraft-comet distance, A, 
was less than 10 4 km. That the HRI detected these point sources is significant. The HRI 
visible camera (HRIVIS) is not in optimal focus, and has an accordingly large and doughnut- 
shaped point spread function (PSF) ( Lindler et al.|2007 ). Thus, unresolved sources are easily 
discriminated from cosmic ray impacts and hot pixels. None of the point sources appear to 
be resolved, therefore they must be smaller than the FWHM of the HRIVIS deconvolved 
PSF (approximately 3 m at A = 700 km). 

In addition to point sources, streaks of various lengths are present in many of the images. 
The lengths of the streaks are correlated with exposure time (longer exposures produce longer 
streaks), and inversely proportional to instrument pixel scale (streaks in HRIVIS images are 
longer than in MRI). The streaks are parallel to each other, and the orientation is consistent 
with the spacecraft's motion as it tracks 103P's nucleus. Stars and other distant objects 



produce streaks of uniform length and orientation, but we find a variety of streak lengths 
in a single image, down to a few pixels long, which indicates they are spatially close to the 
nucleus and not celestial sources. 

Altogether, the above observations indicate that the nucleus of Comet 103P is surrounded 
by a coma of thousands of large particles. We will use the term "particle" to refer to 



the observed point sources in this paper. A'Hearn et al. (2011) demonstrated that these 



particles are centimeter-sized or larger. Large particles are not unique to this comet, but 
these observations are unique to large particles. These are the first images of sub-meter 
particles from a comet seen as individual objects. Prior images in visible/infrared light 
(e.g., comet dust trails; Sykes and Walker 1992, Ishiguro et al. 2002, Reach et al. 2007) 



and with radar (Harmon et al.]|2004) probed collections of many large particles. The only 



other observations of large individual particles from comets that we can think of are those 
of meteors (from meteor streams associated with particular comets) and milligram particle 
impacts on spacecraft. Although the latter examples may only be of order 0.1-1 mm in size. 
In the following paper, we present our methods for detecting and measuring the photo- 
metric properties of the large particles in the coma of Comet 103P. We will convert their 
fluxes into sizes, discuss their spatial distribution, attempt to constrain their composition, 
and compare them to large particles observed in radar observations of the comet. Assuming 
they are icy, we also constrain their contribution to the comet's total water production rate. 



2. Observations 

Images of Comet 103p Were taken with Deep Impact's Medium Resolution Instrument 
and High Resolution Instrument CCD cameras. Both cameras have 1024 x 1024 pixel arrays; 
the HRIVIS pixel scale is 0.413" pixel" 1 (2 u rad), and the MRI pixe l scale is 2.06" pixel" 1 
(10 urad). The instruments are described by Hampton et al. (2005), and their calibration 



by Klaasen et al. (2008, in prep.). The MRI and HRIVIS data are available in the NASA 



Planetary Data System (PDS) archive (McLaughlin et al. 2011a[b ). Images taken with the 
MRI are labeled with the prefix mv, and for the HRIVIS the prefix is hv. We keep this 
convention throughout the paper. We used a pre-PDS data set, but the only significant 
difference is the absolute flux calibration constant, which we account for in our photometry. 
Optimal focus for the HRIVIS instrument occurs about 6 mm before the focal plane. 
For this instrument, it is often beneficial to work with spatially restored (i.e., deconvolved) 
images. We follow the method of Lindler et al. (2007), using the HRIVIS PSF from the 



EPOXI mission ( Barry et al.|2010 ), to restore raw HRIVIS images to near diffraction-limited 
resolution. The HRIVIS images are deconvolved (Fig. [6b) with the Richardson-Lucy (R-L) 
method modified to handle non-Poisson CCD readout noise ( Snyder et al.|199~3 Lindler et al. 
2007). Note that R-L restoration methods conserve flux. 

The large particles are easily seen in MRI and HRIVIS images at A < 10 4 km. Outside 
of this range, the particles become increasingly faint, overwhelmed by the diffuse coma's 
surface brightness, and confused with stars, which are unstreaked point sources when the 
spacecraft-comet distance is large. In a manual search, the earliest image in which we 
have identified particles is HRIVIS image hv5000096, taken 22,400 km from the nucleus 
(CA-30 min, 1500.5 ms exposure time). The first MRI image in which we have positively 
identified particles is mv5002025 (8580 km, CA-11 min, 500.5 ms). 



Most images at A < 10 4 km use the CLEAR1 filter. For sources with solar-like spectra, 
the CLEAR1 filters have effective wavelengths of A e = 0.625 um for HRIVIS, and 0.610 um for 
MRI, and a FWHM of ~ 0.5 um. The absolute calibration uncertainties for the CLEAR1 



filters are 5% for HRIVIS, and 10% for MRI (Klaasen et al.l 2008). Large particles have 



not yet been identified in the HRI infrared spectrometer (HRIIR) scans (Protopapa et al. 



2011), but they may be discovered in the future when we can predict the location of the 



brightest particles in IR data using their MRI and HRIVIS derived 3D positions and velocities 



(Hermalyn et al. this issue). 



A summary of all images used in this paper is presented in Table [TJ Note that the HRIVIS 
images used in this paper are only a small sub-set of the HRIVIS particle observations. We 
have limited our investigation to the closest-approach images because of the great amount 
of time that must be taken to analyze the HRIVIS data. We also do not use any of the 
50.5 ms MRI images in the PDS archive. These images were taken along with the 500.5 ms 
images, but few particles are found in these images due to the decreased sensitivity and great 
distances from the nucleus. 

3. Detection, photometry, and completeness 

3.1. MRI images 

We automatically detect and measure all point sources throughout all MRI images taken 
within A < 5200 km of the comet nucleus using DAOPHOT as included in the IRAF software 



package (Tody 1993). In summary, we: (1) estimate the PSF of the instrument and verify 
that it accurately measures particle fluxes; (2) remove the diffuse coma from the images; 
(3) mask the nucleus and bright jets; (4) additionally mask stars in images at large A; 
(5) measure the completeness of our particle census (i.e., detection efficiency); (6) estimate 
the photometric uncertainties, and correct for any biases; and, (7) clean bad PSF fits and 
crowded regions from our photometry lists. Below, we provide details on our methods. 

Particles were detected by searching for point sources (FWHM = 1.5 pixels) with a 
peak amplitude greater than 5<r above the local background: la per pixel rj 1.1 x 10~ 3 , 
0.52 x 10~ 3 , and 0.13 x 10~ 3 W m~ 2 um -1 sr _1 for 41, 121, and 501 ms frames, respectively. 
Twenty-five isolated and bright point sources were selected by hand from image mv5004056 
to estimate the MRI PSF. The sources were fit by DAOPHOT with a variety of functions, 



and the best-fit PSF used a Moffat function (/3 parameter of 1.5; Moffat 1969) with a look- 
up table of empirical residuals. We compared DAOPHOT-derived fluxes of bright isolated 
particles to those measured with aperture photometry. The values agreed, with a mean 
error of AF\/F\ = —1.9%, where AF\ is the difference between the DAOPHOT flux and 
the aperture flux (F AjDAO phot - -Fx.aper)- 

DAOPHOT iteratively estimates the background as it fits each point source and groups of 
point sources, but we found — especially for detections closer to the nucleus — that subtracting 
an initial estimate of the background improved DAOPHOT's ability to find and measure 
sources. Furthermore, the MRI CCD has an instrumental background that varies row- 
by-row at the 1-4 DN level (the MRI and HRIVIS calibration constants for CLEAR1 are 
3.527 x 10" 5 and 1.209 x 10~ 4 W m" 2 urn" 1 sr" 1 DN" 1 s, respectively). The EPOXI pipeline 
includes a routine to detect and remove this background, but it is not executed on images 
dominated by coma, such as those in our analysis. We derive our background estimate for 



each image using two median niters (i.e., low-pass niters). We first apply a 3 x 25 pixel (row 
x column) filter that defines the background stripes, as well as much of the coma. A second 
25 x 3 pixel filter defines most of the residual coma, especially where it is extended along 
image columns. A median-filter subtracted image is presented in Fig. [T|d. 

Bright jets and the nucleus complicate the PSF fitting. Therefore, we mask these regions 
on the image. The mask is derived from a morphological gradient with an 11 x 11 pixel 
uniform structuring element, which effectively generates an image of peak-to-peak values in 
a moving 11 x 11 pixel box (Fig. lib). We threshold the gradient image at a value of 0.1 
W m -2 urn -1 sr _1 , then dilate the mask by 30 pixels to fill small gaps between features, and 
cover some of the nearby jets (Fig. [Iji). 

DAOPHOT will occasionally fit streaks in the MRI images with PSFs. Most of these 
fits have poor residuals, and are removed before the final analysis (described below). But 
we found several streaks persisting into the final analysis in images taken at large distances 
from the comet (A > 5000 km). These streaks were stars and not particles. DAOPHOT 
seemed more likely to fit streaks in these images because the stars are much brighter than 
the particles, and are only streaked several pixels or less. We generated an additional mask 
by searching for elongated sets of pixels above the background in order to remove streaked 
celestial objects from these images before processing with DAOPHOT. At closer comet- 
spacecraft distances, stars are streaked over a greater number of pixels (up to a few hundred 
pixels at closest approach), and are easily rejected without the mask. 

We verified the completeness (i.e., detection efficiency) of our method by inserting artifi- 
cial point sources into each image, re-executing our scripts on the new images, and examining 
the output to determine if the artificial sources were detected. The completeness is then the 
fraction of detected artificial sources per flux bin. We inserted 1% of the number of detected 
point sources (with a minimum of 10 new particles) uniformly over the image but avoiding 
the nucleus. The fluxes were picked from a distribution uniform in log-space, based on the 
minimum and maximum fluxes measured in the image. The process is repeated, each time 
starting with the original image, until 3000 artificial particles have been tested. Examples 
of our MRI completeness test results are shown in Fig. [2} The completeness is directly cor- 
related with exposure time; in general, it is easier to find and measure fainter particles in 
images with longer exposures. In our MRI population studies, we will only consider particles 
brighter than the 80% completeness level for each image (listed in online Table [2] as i^mm)- 

The completeness tests can be used to derive the photometric uncertainties. As an 
example, we plot the input flux versus the DAOPHOT flux error for all artificial particles 
added to image mv5004046 (Figwl). In this image, the mean flux error AF\/ F\ = —8%, 
and the mean error for all images ranges from —10% to —2%. However, the distribution is 
not Gaussian as there is a significant tail at large negative errors. If we instead take the 
resistant mean (clipping at 3er), we find the mean AF\/F\ ranges from —3.5% to —0.7%, 
with per image standard deviations ranging from 6% to 10%. As discussed above, we also 
found a slight negative error when we compared the DAOPHOT photometry to aperture 
photometry of isolated bright sources. To account for this effect, we will increase the fluxes 
of all particles by 2%. In addition, we adopt 8% as our particle flux uncertainty. 

The large flux errors (\AF/F\ > 10%) are likely misidentified artificial particles, artificial 
particles that have merged with brighter particles in the PSF fitting process, and PSF 
over /under fitting. Any particles with significantly large flux errors must be identified and 
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Figure 1: The process by which we remove the background and mask the nucleus in MRI images, (a) A 
fully calibrated image, mv5004046. (b) The median-filter subtracted image, (c) Image of the morphological 
gradient of image b. (d) The final mask, defined where image c > 0.1 W m~ 2 urn -1 sr . The mask has 
been dilated by 30 pixels to fill gaps and cover some additional coma. The image scale and orientation are 
given in panel d. The projected vectors are comet-Sun (©), comet-Earth (©), and spacecraft velocity (i>di)- 
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Figure 2: Example point source detection completeness functions. All images were taken within A = 
694 - 704 km, except mv5004009 (501 ms) taken at A = 1858 km. The mv5004046 sub-frame is chosen to 
match the same area as covered in hv5004031. The HRIVIS test particles are selected from one of 9 specific 
fluxes (marked by squares and circles), whereas the MRI test particles are selected at random from a range 
of fluxes. The lowest of the HRIVIS bins is 10,000 e". 



removed from our photometry lists so that they do not affect the final results. We clean the 
DAOPHOT output by rejecting: (1) bad PSF fits as identified by DAOPHOT (i.e., reduced 
X 2 > 2.0); and (2) any particles with a total residual in a 3 x 3 pixel box greater than 3cr 
from the background. For the 41 ms images 62 — 78% and — 2% of the detected sources 
were rejected by the two criteria, respectively. For the 121 ms images the rejection rates were 
5-32% and 2-5%; for 501 ms, the rates were 0-1%, and 3-20%. Occasionally DAOPHOT 
fits portions of the particle streaks, and these criteria help remove such spurious fits from 
the photometry lists. Examples of point source cleaned images are shown in Figs. [4] and [5j 
After point source cleaning, sources not fit by DAOPHOT are evident. These particles can 
be accounted for in a statistical sense with the image's completeness function, e.g., when we 
are displaying the flux distribution of the particles, we divide the measured flux distribution 
by the completeness curve to show the true flux distribution. 

In Fig [5] we specifically show a region of high point source density, and the large number of 
bad-PSF fits resulting from point source crowding near the nucleus (> 0.025 pixel -1 ). These 
sources are primarily rejected by the y 2 criterion. Regions with a rejection rate > 20% in a 
25 x 25 pixel area will be masked from our analyses. 

3.2. HRIVIS images 

The analysis of the HRIVIS images required an approach independent from the MRI 
analysis. Figure |6^l shows an HRIVIS image of the large particles and Fig. |6Jd is the restored 
(i.e., deconvolved) image. The broad defocused PSF and the longer trails resulting from 
parallax motion make PSF fitting even more difficult than in the MRI images. Our approach 
is to instead detect and measure a particle's brightness in restored images. In summary, we: 
(1) iteratively remove the background from each image; (2) detect point sources and streaks 
in the restored images; (3) measure particle fluxes; and, (4) measure the completeness of our 
particle detection scheme. Below we detail our methods. 

We use an iterative technique to remove the background before final object detection 
and photometry. The initial background is set to the deconvolved image filtered with a 
51 x 51 pixel median filter. We then fit a smooth surface (using cubic splines) to the filtered 
image to obtain our background image. The objects within the image will bias the median 
filter and the resulting image overestimates the background. To minimize this bias, we repeat 
the process after ignoring all pixels more than 400 e~ above the background (object pixels) 



during median filtering (1 e~ s _1 = 4.32 x 10~ 6 W m -2 urn -1 sr; Klaasen et al. in prep.). 
We repeat this process three times to obtain and remove our final background (Fig. |6b). 

We detect the objects in the image by creating a mask of all pixels more than 1,000 e~ 
above the background level (Fig.[6Ji). We selected the 1,000 e~ threshold by visual inspection 
of the restored image to avoid detection of restoration artifacts that result from "ringing' 



)> 



around the brighter objects (Lindler et al. this issue). We define contiguous masked pixels 



as part of the same object. Visual inspection of this mask shows many instances where a 
trailed particle near the detection threshold is split into multiple objects because of noise 
and potential variations in particle brightness during the time of the exposure. To avoid this 
problem, we examine each detected object (starting with the brightest) and adding all pixels 
within a factor of 3 of the brightest pixel in the object. We define the resulting contiguous 
dark regions as our detected objects (Fig. [6l 
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Figure 3: Input flux, F\, plotted versus DAOPHOT error, AF\/F\, for all artificial particles inserted in 
image mv5004046. Most fluxes are measured to within 8% (la) of their true value, with a slight skew toward 
negative errors. 
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Figure 4: (a) A detail of image mv5004046, after background coma subtraction, (b) The same image after 
being cleaned of point sources detected and fit by DAOPHOT. In both images, particles with large residuals 
or low signal-to- noise flux ratios have been circled. Similar to Figure 111 the image scale and orientation are 
shown. The square in the lower right image shows the field of view of panels a and b with respect to the 
nucleus. 
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Figure 5: The same as Fig. [4j but for a region of high point source density in image mv5004029. 
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Figure 6: (a) Original HRIVIS image hv5004030. (b) Image deconvolved with the R-L method, (c) Image b 
after background subtraction, (d) Mask of pixels greater than 1000 e~ above the background, (e) Mask of 
detected sources, (f ) Mask of the sources used to compute the flux distribution slope (here, length 5 pixels 
or smaller and total flux greater than 20,000 c~). 
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We perform photometry on the detected particles by starting with the particle mask, and 
adding the flux from their nearest-neighbors (i.e., we grow the mask by 1 pixel, and sum 
the pixels together). Non-linearities are a concern when restoring an image with non-linear 



deconvolution algorithms (Lindler et al. 1994). Fainter point sources will be broader than 



brighter point sources in a restored image. We can minimize any non-linearities by increasing 
the size of the region used for photometry. However, this results in more contamination of our 
results by nearby sources. To test for non-linearity we increased the size of the photometry 
regions by adding additional neighbors (i.e., including next-nearest neighbors, etc.). These 
larger regions showed that non-linearity artifacts are absent from particles brighter than 
10,000 e" (F\ = 1.4 x 10~ 12 and 4.6 x 10~ 13 W m~ 2 urn -1 for exposure times of 126 and 
376 ms). On average, the effective aperture for point sources is approximately a 3 pixel 
radius circle. 

As done with the MRI analysis, we measure the completeness by inserting artificial 
objects in the images of varying brightness and parallax length. These artificial objects are 
added to the raw images before deconvolution. The completeness test also gives additional 
evidence that our measured fluxes are linear with object brightness in the range of brightness 
being considered. For the longer trails, the completeness decreased rapidly with flux. We 
therefore decided to ignore particles with a length more than 5 pixels in the restored image. 
In effect, this criterion limits the line-of-sight distance of the particles from the nucleus to 
< 1 — 2 km at closest approach (A = 694 km), and < 2 — 3 km at A = 900 km. If the flux 
distribution is independent with distance from the nucleus, then the relative distribution 
of number versus flux will remain unchanged. Limiting the length of the objects also has 
added benefits. Long streaks are much more likely to be a combination of multiple particles 
and their photometry is also more strongly affected by errors in the background estimation 
(they cover more pixels). Examples of our HRIVIS completeness test results are presented 
in Figure [2] There is little difference between the 126 ms and 376 ms completeness curves in 
the example. Therefore our completeness is not as strongly affected by background noise as 
it is by other factors, e.g., point source crowding and deconvolution artifacts. In addition, 
our HRIVIS photometry is generally limited by the linearity criterion of F > 10, 000 e~. 

4. Particle fluxes 

4-1. Flux distribution 

We measure the fluxes of all particles and compute their flux distributions, dn/dF. The 
particle fluxes approximately follow a power-law distribution. Therefore, we fit each flux 



distribution with a function of the form dn/dF ex F a using the method by Clauset et al. 



(2009), which is based on Kolmogorov-Smirnov tests and maximum likelihood fitting. We 
restricted the flux range to particles brighter than our 80% completeness estimate. This 
restriction increases the statistical errors in the fits but reduces possible bias errors resulting 
from photometry contamination from neighboring particles and residual background. More- 
over, the fitting method does not include a correction for completeness, so restricting the 
fluxes helps mitigate the effect of incomplete counting; the maximum completeness correc- 
tion over our fit range is 20%, which is insignificant in comparison to the strong power-law 
slopes (see Fig. S6 of A'Hearn et al.||2011). Our results are listed in Table |2J and summarized 



in Table [3j In Figs. 7\ and |8j we plot the final flux distributions and their best-fit trends. 
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The fluxes are plotted as a function of F\j Q0 = F\A 2 /700 2 so that they may be directly 
compared to each other (variation of the phase angle is small, ranging from 79 to 92°, and is 
therefore not considered). The error- weighted mean slopes for the MRI and HRIVIS frames 
are -3.74 ± 0.02 and -2.85 ± 0.03, respectively. 

The MRI and HRIVIS slopes are significantly different from each other. To verify this 
difference, we also computed the flux distributions from MRI sub-frames of images chosen 
close in time to the HRIVIS images. The relative position of the HRIVIS and MRI fields 
are constant with time, therefore measuring the flux distribution in the same field of view is 
straightforward. The comparison, however, is not exact because the frame-to-frame parallax 
of individual particles can be significant (Hermalyn et al.|this issue). After bad-PSF rejection 



and completeness tests only « 130 particles could be used to determine the flux distributions 
in the MRI sub-frames. This is a factor of 3-10 fewer than the number of particles in the 
HRIVIS images. The difference primarily relies on the larger aperture of the HRI primary, 
which allows for a fainter point source detection limit due to the finer resolution and greater 
light gathering area. Our best-fit parameters for the MRI sub-frames are listed in Tableland 
summarized in Table [3J The uncertainties in the power-law slopes are quite large, ranging 
from 0.18 to 0.31, owing to the paucity of particles available in the MRI sub-frame images. 
The weighted mean slope is —3.59 ± 0.07, which is shallower than the mean MRI full-frame 
slope, but they agree at the 2<r level. 

Taken together, the mean MRI full-frame, MRI sub-frame, and HRIVIS best-fit power- 
law slopes show evidence that the flux distribution steepens with increasing flux. Inspection 
of the MRI full-frame distributions in Fig. [7] and the Kolmogorov-Smirnov probability (Pks) 
values reveals that a power-law size distribution is not a good fit to most of the flux dis- 
tributions. This poor fit is especially true for the 121 ms frames at A < 900 km where 
there is a significant curvature in the log-flux distributions of these images. The curvature 
is non-existent or not prevalent in the shorter exposure MRI data, the distant MRI data 
(A > 1000 km), and the HRIVIS data. 

It is important to recognize two points. First, the HRIVIS distributions are poorly 
populated at -F\,7oo > 10~ n W m~ 2 urn" 1 , whereas the 121 ms MRI distributions are well 
defined up to -Fa,7oo ~ 3 x 10 -11 W m~ 2 um -1 . This difference is simply a matter of the 
larger solid angle observed by the MRI camera as compared to the HRIVIS camera; a larger 
solid angle allows for more of the brightest particles to be detected. Second, the HRIVIS 
distributions are valid to lower fluxes than the MRI distributions (column i^mm in Table [3]), 
because the larger primary of the HRI allows for fainter particles to be accurately measured. 
Taken together, these differences demonstrate that the two instruments are measuring two 
different particle flux regimes. In Fig. [9| we plot our best-fit power-law slopes versus the 
minimum flux used in each fit. A trend is evident; steeper slopes are correlated with larger 
-^A,min, but due to the different techniques, fit uncertainties, and fields of view involved, 
the correlation is not conclusive. We remind the reader that the correlation is not due to 
incomplete photometry because we restricted our data sets to fluxes where the completeness 
is better than 80%. Given this restriction, none of our completeness corrections are strong 
enough to have a significant effect on the best-fit slopes. 

In the absence of any additional information on the flux distribution we will proceed 
assuming that the differences between the MRI and HRIVIS flux distributions reflect the true 
flux distribution of the entire particle population. As alternatives to a power-law distribution, 
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Figure 7: Flux distributions (dn/d\ogF) and their best- fit power-law functions for all full- frame MRI data 
sets listed in Table [2j The flux distributions have been corrected for completeness, and their error bars are 
based on Poisson statistics. The flux distributions are grouped by their integration times, and sorted by 
time with the pre-closest approach images at the bottom. Each i-th flux distribution (i = 0, 1, . . .) has been 
offset along the y-axis by 10* , and labeled with A for clarity. Additionally, the fluxes have been scaled by 
(A/700) 2 so that they are directly comparable to each other. The gray x symbols mark particles not used 
in the fit. 
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Figure 8: Same as Fig.[7J but for all HRIVIS images listed in Table [2) 
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Figure 9: Best-fit power-law slope versus minimum flux used in the fit for all images listed in Table [2j 
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we also considered a power-law with an exponential cut off, and a double power-law. The 
functional forms are 

din 

— oc F* exp (-F/71) , (1) 

dn F h ( F\ M2 
rfF 72 V 72/ 

where /% are power-law slopes, and 71 is the flux at which the exponential decay reaches 
a factor of exp(— 1) = 0.37 , and 72 is the turn-over flux from the low- to the high-flux 
slopes (both 7 parameters are specified for A = 700 km). The double power-law fits did not 
converge on a single solution. The power-law with exponential cut-off function, however, is 
a better fit to the closest-approach data (A < 1000 km), with a mean slope and cut-off equal 
to -3.18 ± 0.07 and (5.7 ± 0.3) x 10 -12 W m -2 urn -1 . Unfortunately, these best-fit values 
are poor fits at larger distances where the flux distribution better agrees with a power-law 



(Fig. 10). We suggest that the 121 ms distributions have a systematic effect that causes 
under- counting of the largest flux bins. Again, we note that our completeness corrections 
are not strong enough to significantly affect our fits, and that the distributions as shown in 



Figs. [7] and 10 have been corrected for completeness, yet the curvature remains in the closest 



approach data. 

To test the hypothesis that the power-law slope is affected by the source density, we 
generated completely synthetic data sets with DAOPHOT, based on our MRI PSF. We 
created 8 images each with the same power-law flux distribution (dn/dF oc F -3 ' 5 ) and 
maximum point source flux (FA )max = 1 x 10 -10 W m -2 urn -1 ), but F\ >m i n varied from 
5x 10 -12 down to 1.5 x 10 -13 W m -2 urn -1 . We inserted the same number of point sources per 
logarithmic flux bin, but because each image's flux distribution spanned a different range, the 
total point source density varied from 4 x 10 -4 pixel -1 up to 2.6 pixel -1 . We then measured 
each image's flux distribution with DAOPHOT in a manner similar to our MRI images. 
Once the input column density reached 0.08 pixel -1 , our best-fit slopes steepened from the 
input —3.5 down to a minimum of —4.7 at the highest densities. Thus, the rollover in the 
flux distribution seen in the 121 ms MRI images at A < 1200 km could be a manifestation 
of this effect. 

The slopes derived from the distant MRI and the 41 ms MRI images are consistent: 
a = —3.80 ± 0.02 and —3.83 ± 0.05, respectively. They agree despite the wide range in the 
number of particles fit (100-1000) and spacecraft-comet distances (900-5000 km, although 
the best constraints are within 3000 km). Therefore, we consider these results to be robust 
and representative of the true flux distribution for i*\,7oo ^ 1 x 10 -11 W m -2 um -1 . The 
shallower HRIVIS slopes appear valid for 1 x 10 -12 < Fx < 1 x 10 -11 . We see no reason to 
trust one data set over the other, or to assume that a single power-law would be valid for all 
fluxes measured. Therefore, we will use a broken power-law for the remainder of the paper: 



dn ,_,_, 


a = 


-2.85 


-J^^F 


) 




dF 


1 a = 


-3.80 



for F Ai700 < 7.2 x 10 -12 W m -2 um -1 

(3) 



A,700 

for F Ai70 o > 7.2 x 10 -12 W m -2 um -1 



where -Fa, 700 = 7.2 x 10 12 is the break in the power-law, —2.85 ±0.02 is the average HRIVIS 
slope, and —3.80 ± 0.02 is the average MRI slope, measured from images at A > 900 km 
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Figure 10: Best-fit power-law with exponential cut off, derived from 121 ms MRI images taken within 
A = 1000 km, compared to selected flux distributions, labeled with "A (exposure time)" in km and ms. The 
histograms have been offset by a factor of 10* for clarity. 
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Figure 11: Combined HRIVIS and MRI flux distributions (F\ > F A ,min), scaled to N(F X = 7.2 x 10~ 12 ) = 1. 
For the MRI data, only full-frame distributions from A > 900 km are shown. The solid line is our best-fit 
broken power-law (Eq. p5j). 



to avoid the apparent crowding effects at closest approach. The break was derived by least- 
squares fitting the combined MRI and HRIVIS data sets. The broken power-law and the 



combined HRIVIS and MRI flux distributions are presented in Fig. 11 This model is a good 
match to the ensemble (HRIVIS + MRI) data set. 

4-2. Total particle flux 

Due to the great numbers of large particles, we can only accurately measure the total 
scattered flux from the brightest particles in each image. Our best estimate of the frac- 
tion of the coma flux attributable to large particles is 0.11% for fluxes ranging 0.9 — 4.7 x 
10~ n W cm' 2 pin -1 (row labeled "MRI (distant)" in Table pi). This estimate was derived 
from images with relatively little point source crowding (A > 2000 km), and a moderate 
number of particle measurements (N > 300 particles). In order to build a complete cen- 
sus of the particles, we need to extrapolate the results from that limited range down to 
the faint end of the distribution. We find that the brightest particles observed throughout 
the closest approach images are consistently near -Fa, 700 — 4.5 x 10 -11 W cm -2 pin -1 . On 
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the other end, the faintest particles we can manually find and measure have fluxes of order 



^\,7oo ~ 10 13 W cm 2 um l (Fig. 12). If the very steep flux distribution we derived in £4.1 



extends from -F\,7oo = 1 x 1CT 12 down to 1 x 1CT 13 W cm -2 um -1 , there should be several 
million particles per image. However, at such great pixel densities (1 particle per 17 pixel 
area) it should be very difficult to find faint isolated particles (the core of the HRIVIS native 
PSF has an area of 113 pixels). Furthermore, if we let the flux distribution continue down 
to 1 x 10~ 14 W cm -2 um -1 , the large particles account for 100% of the coma flux, leaving no 
room for any fainter particles (millimeter sized and smaller). Therefore, the flux distribution 
must change or be truncated at fluxes fainter than .F\,7oo = 1 x 10~ 12 W cm~ 2 um -1 . The 
uncertainty in the lower flux limit is a major source of error in all estimates of the total 
number, flux, cross section, and similar derived quantities. We list the results from our ex- 
trapolations in Table |ij We estimate that the large particles account for 2-14% of the total 
coma flux near the nucleus. 

If we take the the first image in which particles are easily seen, mv5002051 at A = 5148 km 



(Fig. 13), we can estimate the total number of particles within a 20.6 km radius (the largest 
circular aperture centered on the nucleus that fits within the image). We find a total coma 
flux of (2.15 ± 0.22) x 10 -7 W m~ 2 um -1 (includes the bright jets, but not the nucleus) and 
assume that 2-14% of that flux is attributable to large particles. The results are presented 
in Table [4] and will be used below when we estimate the cross section, mass, and water 
production rate of the large particles. 

Up to now, we have not addressed streaked particles in our MRI images. If a significant 
fraction of the flux in our images at A w 2400 km is contained within streaked particles 
then our estimates on the total flux and number of large particles will be systematically low. 
However, we find that few particles are streaked and the correction to include any possible 
streaks is negligible. To demonstrate, we developed a Monte Carlo simulation that uses the 
spacecraft position and velocity to estimate the volume of space that is smeared. Within a 
20.6 km projected radius, 2% the field-of-view is smeared over > 1.5 MRI pixels, and only 
particles outside of 120 km are smeared this much. If the particle density follows a r~ 2 
profile, the fraction of smeared particles reduces to effectively zero. 

5. Particle size and composition 

In the absence of any compositional information on these large particles, we assume two 
cases to demonstrate their likely range of sizes. First, we will consider that the particles are 
refractory, and photometrically behave like comet nuclei (the "dusty" case). Then, we will 
consider that the particles are icy, and behave like the icy satellite Europa (the "icy" case). 
We stress that these two cases are examples only. They may not reflect the true nature 
of the particles, but they do yield useful limits on the particle sizes. We will show that 
the particles are much larger than the 0.1-100 um sizes typically considered in comet dust. 
We purposefully avoid interpreting the large particles with phase functions that have been 
derived for comet comae. Light scattered by comet dust comae is expected to be dominated 



by dust grains in the sub-micrometer to micrometer size range (Kolokolova et al. 2004), 
which scatter optical light differently than centimeter-sized particles. 

First we take the case in which the particles photometrically behave like comet nuclei, 
i.e., they have a very low albedo, and have a phase angle behavior like a macroscopic object 
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Figure 12: Two sub-frames of HRIVIS image hv5004031 (A = 694 km) showing the crowded field of particles: 
a) an MM context image (mv5004046) with the approximate HRIVIS field of view outlined with a box (image 
is log scaled from 0.001 to 1.0 W m~ 2 urn -1 sr _1 ); b) HRIVIS image, sub-frames c/d and e/f are outlined and 
labeled; c) HRIVIS sub-frame; d) the same as c, but deconvolved to enhance the spatial resolution; e) another 
HRIVIS sub-frame; f) panel e, deconvolved. Two particles with fluxes near 1 — 2 x 10~ 13 W m~ 2 urn -1 have 
been circled. The HRIVIS field of view is 1.4 km, and the sub- frames are 0.42 km. 
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Figure 13: Image mv5002051, background coma removed. Point sources identified by DAOPHOT with fluxes 
greater larger than the image's 80% completeness limit are circled to show the extent of the large particle 
coma (particles below this limit can still be seen by eye). An image of the nucleus has been inserted into 
the central masked region. The box marks the location of the inset image. Stars, most of which have been 
automatically masked, are streaked over 13 pixels in this image. 
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and not like small dust grains. We refer to this case as the "dusty case," but, just like 
comet nuclei, these model particles may be internally icy. We adopt the geometric albedo 
(A p = 0.049) and phase function (— 2. 5 log $ = 0.0460 mag for 9 in units of degrees) of 
Hartley 2's nucleus (Li et al. this issue). The geometric albedo is defined as the ratio of 



the energy scattered from the object toward a phase angle of 0° to that scattered from a 
white Lambertian disk with the same cross section (cf. Hanner et al.|1981 for this and other 
relevant albedo definitions). 

1 x 10~ 13 



In Section 4.1, we found that individual particle fluxes range from F; 



A,700 = 



to 



4.5 x 10 n W m 2 pm 1 , where F\j 00 is the particle flux normalized to a distance of 700 km. 
To convert between flux and cross section, we use the formula 



F x = A p $(9)a 



S 



A,© 



A 2 



(4) 



where F\ is the particle flux in units of W m~ 2 um -1 , a is the cross-sectional area of the 
particle (cm 2 ), S\ tQ is the solar flux density at 1 AU (W m -2 um -1 ), r^ is the heliocentric 
distance of the particle (1.064 AU), and A is the spacecraft-particle distance (cm). For the 
HRIVIS and MRI CLEAR1 filters, we use a solar flux density of 1471 and 1435 W m" 2 pm" 1 , 
respectively (Klaasen et al. in prep.). To convert from cross section to effective radius, 
we assume a spherical geometry: a = iia 2 . Altogether, assuming a comet nucleus-like 
photometric behavior, the effective radii of the particles range from 10 to 221 cm. The largest 
of these particles (4 m diameter) are just over the resolution of the HRIVIS reconstructed 
frames (about 3 m at a distance of 700 m). We consider these sizes to be an upper limit to 
the true particle sizes. For a lower-limit estimate, we consider the case where the particle 
scattering function is similar to the icy satellite Europa: A p = 0.67 and $ = 1.0 — 0.010 + 



2.2 x 10~ b 9 2 (Buratti and Veverka 1983, Grundy et al. 2007). At a phase angle of 80° 



Europa's phase function is 0.34. For this parameter set the effective radii are 0.8 to 17 cm. 

Following Eq. [IJ where flux is proportional to cross-sectional area, we can compute the 
total observed particle cross section and add it to Table |4j We have assumed our icy particle 
case, using the photometric parameters of Europa. To instead use our dusty case, multiply 
the cross sections in the table by 158. 

Our flux distributions imply a very steep size distribution. Assuming a spherical geome- 
try, the differential size distribution is 



dn 
da 



2N F? +1 a 2a+1 , 



(5) 



where F\ is the flux from a 1 cm radius particle in units of W m 2 pm 1 , and particle radius 
a is in units of cm. The constant N is the solution to the equation: 



N = N, 



*A,. 



" A, min 



dn 



dF, 



(6) 



using the appropriate values from Table |4j Our low-flux power-law slope (a = —2.85) 
corresponds to a size distribution proportional to a -4,7 . This slope is steeper than the many 
size distribution estimates of small through large dust grains (a ~ 1 pm to 1 mm with slopes 
near —4 to —3) based on grain thermal emission (Lisse et al.|1998, Harker et al.|2002, 2011), 
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grain dynamics QFulle|2004[ [Reach et al.|2007[|Kelley et al.|2008[|Vaubaillon and Reach|20lo| , 
and dust flux monitors on spacecraft ( McDonnell et al.|1987 , Green et al.|2004 ). Specifically 
for Hartley 2, Bauer et al. (2011) estimate the comet's overall size distribution to follow a 



power-law slope of —4.0 ±0.3, derived by comparing .R-band and WISE 12 and 22 um fluxes. 



Epifani et al. (2001) fit an ISOCAM 15 um image of the comet taken about 10 days after 



perihelion with a dust dynamical model. They report a time-averaged power-law slope of 
—3.2 ± 0.1, but inspection of their Fig. 10 suggests —3.8 is more appropriate (their best-fit 
slopes never fall above —3.5). These slopes are shallower than our value of —4.7, but there 
is no requirement that they be the same. 

To better understand the size, and thereby the composition, of the large particles, we 
investigate the largest particle that may be lifted from the nucleus, a crit . This parameter is 
estimated by comparing the force of gravity to the gas drag force at the surface of the comet. 



Meech and Svoreh (2004) integrated the equation of motion for spherical particles ejected 



from a spherical nucleus and found 



^cr it 



9pm n Qv th 
6in 2 ppp N R%G : 



(7) 



where /i is the atomic weight of the gas (amu), m H is the mass of hydrogen (g), Q is the gas 
production rate (s -1 ), v t h is the mean thermal expansion speed of the gas (cm s -1 ), p p is the 
density of the particle (g cm -3 ), p^ is the density of the nucleus (g cm -3 ), Rn is the radius of 
the nucleus (cm), and G is the gravitational constant (cm 3 g -1 s -2 ). With the shape model 



of Comet Hartley 2, Thomas et al. (this issue) estimate the surface gravity of the nucleus to 
be a^v = 0.0019 — 0.0044 cm s -2 , which includes the rotation state of the nucleus. Therefore, 



we re- write the equation from Meech and Svoreh (2004) to use the gravitational acceleration 
at the nucleus 



dcrit 



3pm R Qv th 
lQnp p a N R 2 N 



The bulk material density for dust is ~ 3 g cm -3 and for ice is 1.0 g cm -3 . For the dusty 
case, we will consider porous aggregates of dust with a total density of 0.3 g cm -3 (i.e., 90% 
vacuum). For the icy case, we will at first assume 1.0 g cm -3 , but later consider porous 
aggregates with p p = 0.1 g cm -3 . 

The total water production rate of Comet Hartley 2 near closest approach has been 
estimated via several methods to be Q(H 2 0) « 0.7 — 1.2 x 10 28 s -1 (|A'Hearn et al. 



2011 



Combi et aLl[20TT| [DeUo Russo et al.||20TT| |Meech et al.||20TT| |Mumma et al.||2011j |Knight 



and Schleicher this issue). Yet, the coma contains a significant amount of water ice (A'Hearn 



et al. 2011, Protopapa et al. 2011), which could be supplying a large fraction of the water 



vapor around the comet. Moreover, the water production rate is not uniformly distributed 
over the surface ( A'Hearn et al.|2011 ). We can account for these observations by multiplying 
the water production rate by the ratio f surface/ f active, where f sur face is the fraction of the 
water vapor produced at the surface of the nucleus, and f active is the areal fraction of surface 
that is active. For illustrative purposes, we will assume that 1% of the water vapor is 
produced from 10% of the surface. The remaining 99% of the water vapor sublimates from 
the icy grain halo. 

For a nucleus water production rate of Q (R2O) f surface/ f active = 10 27 s -1 , v t h = 0.5 km s -1 , 
Hartley 2's mean radius (0.58 km), mean surface gravity, and assuming a particle density of 
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pure ice (p p — 1 g cm 3 ), we find a cr u = 8 cm. If we instead take CO2 as the driving gas, 



with a production rate of 10 27 s _1 ( A'Hearn et al.|[20lT ), f surface = 100%, and f active = 10%, 
then a crit becomes 200 cm for solid ice spheres. Based on this exercise, the icy model size 
estimates, a < 17 cm, are reasonable. 

If instead of icy particles, we assume nucleus-like particles with a density of 0.3 g cm -3 , 
our a cr it estimates increase to 28 cm (H 2 0) and 670 cm (CO2). Compared to our size 
estimates of a < 210 cm, dark, dusty particles are plausible if CO2 is the driving gas; water 
can be made consistent if f surface/ f active is increased to 1.0. 

Aggregate particles are more easily lifted by gas drag due to their larger surface area 
per mass. Nakamura and Hidaka (1998) found that the drag force for an aggregate is 



approximately the same as the drag force on an area-equivalent sphere (with an error of less 
than 40% in the large aggregate limit). Since our particle radii are based on the observed 
flux, which is proportional to the cross-sectional area, our radii are already defined by area- 
equivalent spheres. Therefore, by revising our particle density we can use Eq. [8] to estimate 
a crit for aggregates. 

We assembled model particles using a ballistic particle-cluster aggregate (BPCA) method 



(Meakin 1984). As the size of the aggregate grows beyond a few thousand monomers the 



density asymptotically approaches 10% of the bulk material density. Thus, for a refractory 

material with a bulk density near 3 g cm -3 , a centimeter-sized BPCA particle would have a 

density of 0.3 g cm -3 . The a cr u estimates will be the same as in our nucleus-like case above. 

Treating the large particles as aggregates rather than solid spheres better agrees with the 



HRIIR spectra of the coma. A'Hearn et al. (2011) and Protopapa et al. (2011) studied the 



water ice absorption features and found that they are most consistent with icy aggregates 
with monomer radii ~ 1 um. An icy BPCA particle would have a density of 0.1 g cm~ 3 . The 
lowered density for icy aggregates increases a crit by a factor of 10, giving us a healthy margin 
for launching large icy particles off the surface of the nucleus, even where water sublimation 
is driving the activity. 

In summary, the dusty particle case produces very large particle estimates (up to 2 m 
in radius) that are just at the resolution limit of the HRIVIS instrument. Gas expansion 
from CO2 is sufficient to lift these large particles from the surface of the comet if they have 
a comet-like density of ~ 0.3 g cm -3 . The water-ice case produces particle estimates up to 
« 20 cm in radius, which are easily lifted from the nucleus by water or CO2 expansion. 



6. Spatial distribution and origin 

Mapping the spatial distribution of the particles gives clues to their origins and dynamics. 
In Figs. [T4)-[T7|, we plot the column density of particles and total coma surface brightness 



contours for all 501 and 121 ms MRI images listed in Table [2j The column density images 
were derived from our final photometry lists for F\ > F\ jm { n , binned onto a 38 x 38 pixel 
grid. Inspection of the figures reveals that the coma and the particles have different spatial 
distributions. The particles are biased to the anti-sunward direction on scales > 2 — 4 km, 
whereas the coma distribution is dominated by the jets. This asymmetry is especially ap- 
parent in the strong sunward jets, where the particle density is lower than in other regions 
of similar surface brightness. It is not an observational bias; we have masked those regions 
close to the nucleus where the column density and jet morphology interfere with the PSF 
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fitting process. Instead, the low particle density in the sunward jets can be accounted for 
by particle dynamics. We consider three dynamical processes that could be affecting the 
particle distribution: (1) the rotation of the nucleus; (2) solar radiation pressure; and (3) 
a rocket effect from sublimating ice. Hydrodynamic flows from the strong gas production 
rates and asymmetry in outgassing may also play a role in the particle dynamics, but this 
analysis is outside the scope of this paper. 

6.1. Radial expansion and nucleus rotation state 

Rotation of the nucleus has a direct impact on the spatial distribution of particles. To 
estimate this effect, we must first recognize that the particle outflow speeds are low. In 
order for the particles to be seen as point sources in HRIVIS images at closest approach, 

their speeds must be lower than ~ 3 m/0.376 s = 8 m s _1 . A lower constraint is computed by 



Hermalyn et al. (this issue) based on the 3D positions of the particles. They find 0.5-2 m s _1 
to be more typical (note that these speeds are not necessarily radial). At such low speeds, 
the large particles take > 10 3 s to reach 2 km from the nucleus. In contrast, small dust 
grains move much more quickly with outflow speeds expected to be of order 100 m s _1 . The 
fine dust reaches 2 km in as little as 20 s. 

The positions of the major jets are governed by the rotation of the long axis about the 
angular momentum vector, with a period of 18.4 h near closest approach; the long axis is 



inclined to the angular momentum vector by 81° (Belton et al. this issue). Taking 1ms 



as the particle outflow speed from the surface, the nucleus will have rotated 10° by the time 
the particles have traveled 2 km. So, the rotation state has a minor consequence on the 
distribution of particles at 2 km, and their spatial distribution should be closely related to 
their source regions. This observation and the assumption of radial motion suggests that 
the strong sunward jets are not the primary source of the large particles, but instead they 
are ejected from along the long-axis of the nucleus. The jets pointed towards the bottom 



of the images in Fig. 15 would be the next likely source region for particles. There is also 



a large population of particles towards the top of Figs. 14 and 15, but they do not have an 
apparent source region (i.e., this side of the nucleus does not appear to be as active as the 
other regions). However, we note that not all potential sources are apparent in the MRI 



images. For example, the apparent water jet seen in Fig. 5 of A'Hearn et al. (2011), which 
points to the top right of Figs. |16] and 17 has no clear optical counterpart but may also 
contribute to the large particle production (n.b., this water jet does not appear to contain 
ice, and therefore is unlikely to be a source of large icy particles). 

Taking a lower ejection speed of 1 cm s _1 , the nucleus can rotate three times before 
particles travel 2 km. Thus on this length scale, and in the absence of any other perturbing 
forces, the particles would have a spatial distribution correlated with the activity profile of 
the nucleus. Because the CO2 and H 2 gas production rates and the optical light curve peak 
when the small end is pointed towards the Sun ( A'Hearn et al.pOll ), in the absence of other 
forces the large particle density should peak in the solar direction, which is not observed. 

6. 2. Solar radiation pressure 

In principle, solar radiation pressure could redistribute large particles into the anti- 
sunward direction. The acceleration from radiation, a ra d, in units of cm s~ 2 is (Burns et al. 
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Figure 14: Contours of the total coma surface brightness superimposed over particle column density for 501 
and 121 ms MM images mv5002051 through mv5004025. The contours are spaced at factor of 2 intervals, 
and the dashed line is 0.008 W m~ 2 urn -1 sr _1 . The coma image was smoothed with a Gaussian kernel 
function (7 pixel FWHM) before the contours were created. Only particles with F\ > -Fa, mm are considered, 
and the particle bins are 38 x 38 pixels in size. The region closest to the nucleus is masked from the analysis, 
and has been replaced with each epoch's image of the nucleus and inner-most coma. In all images, the 
projected velocity of the spacecraft is toward the bottom, and the sunward direction is approximately to the 
right. 
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Figure 15: Same as Fig. [141 but for images mv5004029 through mv5004053. 
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Figure 16: Same as Fig. [141 but for images mv5004056 through mv5006020. 
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Figure 17: Same as Fig. [141 but for images mv5006024 through mv5006046. 
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where Q pr is the radiation pressure efficiency factor, S Q is the integrated flux density (1.361 x 
10 6 erg cm" 2 s" 1 at 1 AU), a is the geometrical cross section of the particle (cm 2 ), c is the 
speed of light (cm s _1 ), m is the mass of the particle in question, and r^ is the heliocentric 



distance (AU). We have already dropped all velocity dependent terms from a rad (cf. Burns 



et al. 1979). It is common to express radiation pressure with the parameter /3 defined as the 



ratio of the force of solar radiation pressure to the gravitational force from the Sun, 



P 



F rad 5.7 x l(T 5 Q pr 



grav 



p p a 



(10) 



The radiation pressure efficiency is 1 for perfectly absorbing, isotropically emitting spheres. 
For our particles, we again take the icy particle case and derive Q pr from the phase func- 
tion and geometric albedo (van de Hulst 1957): Q pr ~ 1 — A B cosa, where A B is the Bond 
albedo, and cos a describes the anisotropy of the scattered light, where a = 180 — 9 is 
the scattering angle. The relationship between the Bond albedo and geometric albedo is 
Ab = A p J* $(a) sinada (Hanner et al. |l981 ). Altogether, we compute Q pr = 1.58. For the 
mass, we again assume two cases for particles with a = 10 cm: (1) solid spheres with the 
density of ice, 1 g cm -3 ; and (2) compact aggregates of ice (1 um radius monomers) with an 
overall particle density of 0.1 g cm -3 . For the latter case, we note that detailed calculations 
will be required to understand how albedo and the anisotropy of scattering are affected by 
the complex aggregate shape, and we reserve this investigation for future work. We compute 
accelerations of 4.8 x 10" 6 and 0.5 x 10~ 6 cm s~ 2 for 10 cm solid and BPCA aggregate ice 
particles, equivalent to (3 = 9 x 10~ 6 , and 1 x 10~ 6 . For comparison, comet dust trails are 
comprised of grains with j3 < 10~ 3 (Sykes and Walker 1992, Reach et al. 2007). 



In Figs. 14-171 particles are found out to the image edges in the sunward direction. 
However, the sunward/anti-sunward asymmetry is clear on 2-4 km length scales. In the rest 
frame of the comet, the relationship between turnaround distance (d), ejection speed (v e j), 



and acceleration from radiation pressure is d 
yields 



ej 



/2/a ra d- Solving for ejection speed v 



<:i 



v 



oCL\cg pr O ( 



pr OQ 



<■:) 



6.0 x 10" 



,dQ 



pr 



;n) 



for v e j in units of cm s 1 
v ej = 4.3-6. l(pa)- 1 / 2 



2cr l h p p a p p a 

d in cm, p p in g cm -3 , and a in cm. For our icy particle case 
cm s _1 (4-6 cm s" 1 for a 10 cm aggregate, and 1-2 cm s _1 for solid ice). 



If particles are truly ejected at these speeds, solar radiation pressure would take > 10 6 pa s to 
accelerate the particles to the > 50 cm s _1 speeds measured by Hermalyn et al. (this issue). 



This result also implies that the dominant velocity component for particles more than a few 
kilometers from the nucleus would be distinctly in the anti-solar direction, yet only a weak 
asymmetry in the velocity is observed (Hermalyn et al. this issue). Therefore, radiation 



pressure does not govern the dynamics for our icy particle cases. The same conclusion is 
reached for the dusty particle case with a nucleus-like density (Q pr = 1.00, p p = 0.3 g cm" 3 ): 



J ej 



6.3 — 8.9a 1 I 2 cm s 1 or 0.6-0.9 cm s l for a 100 cm particle. 
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If we instead assume an ejection speed of 1 m s _1 for a 10 cm particle, their densities 
must be p p = 2 — 4 x 10~ 4 g cm -3 in order to be turned around by d = 2 — 4 km. Thus, if 
the particles are very fluffy aggregates, they may be ejected at larger speeds, and radiation 
pressure will be able to redistribute them into the anti-solar direction, but they could still 
have instantaneous velocities distributed about the anti-solar vector. 

6.3. Rocket effect 

The rocket effect is the acceleration of the particles due to the sublimation of water ice. 
For a spherical geometry with radial outflow, 

a mt „ = 3 " m » Zt "^ = 1.1 x 10-'» A, (12) 



4p p a p p a 



where Z is the sublimation rate of the particle (mo lee cm 2 s x ), \x is the molecular weight 
of the sublimating ice (18 u for water), mn is the mass of hydrogen (g), and fi ce is the ice 



fraction of the particle (all remaining parameters are in cgs units). Following Reach et al. 



(2009), we define the dimensionless rocket effect parameter a as the ratio of the force from 
the sublimation mass- loss, F TO( ± eU to the force of gravity from the Sun, F grav : 

= F ^L = 2.238 x 10 ^ 3 ^Zv th f tc 

F grav 4GM ePp a V ; 

where G is the gravitational constant (cm 3 g _1 s~ 2 ), and M Q is the mass of the Sun (g). The 
rocket effect parameter should have a heliocentric distance dependency, since the product 



Zvth does not necessarily vary as r^ 2 , but Eq. 13 will serve as a good approximation for 
small Arv The a parameter is analogous to the /3 parameter for dust; both parameters 
quantify a force directed away from the Sun in fractions of the solar gravitational force. 

Assuming no losses from scattering or thermal emission, the maximum water sublima- 
tion rate from ice (Z max ) can be computed from the solar flux density, particle absorption 
efficiency (Q a b S ), and the latent heat of sublimation for water ice (L), 

rv S e Q a b s NA 



rlL 



(14) 



where N& is Avogadro's number, and L ranges from 5.0 x 10 11 to 5.1 x 10 11 erg mol 1 for 



temperatures from 100 to 300 K (Murphy and Koop 2005). For large compact aggregates 



(i.e., porous spheres) and solid particles Q a bs ~ 1, but it is larger for fluffy aggregates since 
they have light scattering properties more like a collection of monomers, rather than a single 



solid particle (Kolokolova et al. 2007). This last comment aside, Z max becomes 1.4 x 10 



18 



molec cm 2 s l for T = 300 K. To account for scattering, Z max will scale with 1 — As, 



where As is the Bond albedo. In £6.2 we computed As = 0.84 for our icy particle model. 
Scattering reduces the Z max of the icy model to Z max ^ cy = 2.2 x 10 17 molec cm~ 2 s -1 . For 
our dusty model, A B = 0.013 and Z max4usty g Z max . 



For comparison, consider the model of Beer et al. (2006) for the sublimation rate of icy 



particles. They employ spherical grains, heated by absorption of sunlight and cooled by 
thermal emission and sublimation, and computed the absorption and emission efficiencies 
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with Mie scattering and effective medium theory considering both pure ice grains, and "dirty- 
ice" grains, i.e., ice mixed with a generic absorber (dust). Mixing dust with the ice has a 
significant effect on the ice equilibrium temperature, but the effect is not a strong function 
of the ice-to-dust mass ratio (they tested m ice /md U st = 0.9 and 0.5). Dust should also affect 
the overall albedo of the particles, but they do not report this parameter. In their Figs. 8 
and 9 they present computed grain lifetimes, defined as the time a grain takes to completely 
sublimate. Their lifetimes for r^ = 1.09 AU are the best examples for our scenario. For 
a > 0.01 cm, pure ice grains have lifetimes of r pure = 1 x 10 10 a s for particle radius a 
measured in cm. Dirty-ice grains have a much shorter lifetime: Thirty = 1 x 10 5 a s. Since 
the large particles spend most of their lifetime larger than 0.01 cm, we can transform their 
lifetimes into sublimation rates: 



'■'dirty 



(15) 



The Beer et al. (2006) model sublimation rates for large particles are Z^ Tiy = 4 x 10 17 



molec cm -2 s -1 , and Z pnre = 4 x 10 12 molec cm -2 s -1 . These values are comparable to or 
less than our maximum sublimation rates. 

Taking Z max ^ cy , a 10 cm radius particle is accelerated at a rate a roc k e t = 0.024 cm s -2 
away from the Sun, yielding a rocket parameter a = 0.042, i.e., the rocket effect is 4% the 
force of solar gravity. For Z = Z max , we compute a = 0.27. Unlike radiation pressure, the 
rocket effect is potentially very strong. 

Acceleration from sublimation will distribute the particles in the anti-solar direction. 
Following our method for radiation pressure, we can constrain the particle ejection speed 
with the implied sublimation rate 



Zd 
p p a 



v 2 ej = 2.2 x 10 -18 — . (16) 



Again, adopting 2-4 km as our typical turnaround distance, and taking Z = Z maXticy , we find 
ejection speeds of 310 — 440a -1 ' 2 cm s -1 for our solid ice case (990 — 1400a -1 ' 2 cm s -1 for 
icy aggregates). These ejection speeds are higher than the instantaneous speeds measured 



by Hermalyn et al. (this issue), but the two speeds do not need to agree since one is at /near 
the surface and the other is out in the coma. For Z = Z max , the resulting ejection speeds 
are increased: v e j = 780 — 1100a -1//2 cm s -1 (solid ice) and 2500 — 3500a -1 / 2 cm s -1 (icy 
aggregates). However, in the above analysis we have assumed that only the sunlit hemisphere 
is sublimating. By distributing the sublimation across more of the surface, we can decrease 
the implied ejection speeds to 10-100 cm s -1 , similar to the speed measured in the coma by 



Hermalyn et al. (this issue). Therefore, we conclude that a sublimation rate excess on the 
sunlit hemisphere of order 10 17 molec cm -2 s -1 readily describes the sunward/anti-sunward 
particle asymmetry. Detailed simulations will be needed to fully account for the observed 
distribution of particle velocities. 

7. Mass and water production rate 

Table [4] includes the total particle mass, based on our icy photometric model and a particle 
density of 1 g cm -3 . For the dusty model and p p = 0.3 g cm -3 , the masses are increased by 
a factor of 672. With our preferred flux lower limit of 0.1-1.0 x 10 -12 , the particle masses 
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correspond to 0.03 - O.IOMat (icy), 23 - 70M N (dusty), where > M N = 2.4 x 10 11 g is the mass 
of the nucleus, assuming a 0.3 g cm -3 density (Thomas et al. this issue). It is clear that the 



dusty case is impossible. The only way we can reduce the estimated total mass for these 
dark particles is by decreasing their densities to well below 10~ 3 g cm -3 . Therefore, we do 
not favor the dusty case, but it does remain as a possible interpretation. Given that the total 



mass lost from the comet per orbit is of order 1% of the nucleus (Thomas et al. this issue), 
even the solid ice particles may be too massive. Porous ice particles (e.g., p p = 0.1 g cm -3 ) 
should be considered the most likely case. 

Icy particles will begin sublimating as soon as they are released from the nucleus and 



warmed by insolation. In £6.3, we computed the rocket effect on the particles due to water 
ice sublimation and concluded that a sublimation rate excess of 10 17 molec cm -2 s _1 on the 
particle's sunlit hemisphere can describe the observed sunward/anti-sunward asymmetry 
in particle column density. We also computed the maximum sublimation rate, based on 
energy balance between absorbed solar radiation and sublimation. We apply this latter 
value, Z max , icy to all of the large particles and compute water production rates. The total 
large particle water production rate within a 20.6 km radius aperture is limited to < 0.6 — 
5 x 10 25 molec s _1 , which is < 0.1 — 0.5% of the total water production rate of the comet 
(kIx 10 28 molec s _1 ). We can change the water production rate by assuming a different 
photometric model as Q oc (1 — A B )(A P &)~ 1 . For our nucleus-like case, we find Q(H 2 0) < 
16 — 80% of the total water production rate, but, unless the particles are fluffy aggregates 
with p < 10~ 3 g cm -3 , we rule out this case based on their mass. The icy particle case yields 
our best estimate of the water production rate, Q(H 2 0) < (1 — 5 x 10~ 3 )Qtotai- 

8. Comparison to other observations 



Harmon et al. (2011) observed comet Hartley 2 with Arecibo S-band (A = 12.6 cm) 
radar at the end of October 2010, about a week before Deep Impact's closest-approach. In 
their average Doppler spectrum, they observe a strong grain-coma echo, with a characteristic 
radial velocity dispersion of 4 m s _1 . The velocity distribution is asymmetric, with a range 
of ~ —50 to +13 m s _1 with respect to the nucleus (negative velocities are away from the 
Earth). The coma has a strongly depolarized echo, which indicates the radii of the largest 



particles are well above the Rayleigh limit of X/2ir = 2 cm. Harmon et al. (2011 ) suggest that 
there exists a significant population of particles with decimeter sizes or larger. We explore 
the possibility that the radar observations may be the large particles imaged by Deep Impact. 
The average radar cross section of the coma was 0.89 km 2 . From our MRI observations, 
we derived a total icy particle cross section of a = 4 x 10~ 4 to 3 x 10~ 3 km 2 within 20.6 km 
from the nucleus. Our cross section is more than two orders of magnitude smaller than 
the radar observed cross section. However, their beam is much larger than what the MRI 
can image when individual particles are detectable. We have found particles out to 40 km 
in Fig. [13J but the Arecibo beam size is 32,000 km at the distance of the comet. With 
speeds of order 4 m s _1 , and lifetimes of order 10 5 s (dirty ice), the large icy particle coma 
would extend out to ~ 400 km. This estimate implies our census of the large particles is 
incomplete by about a factor of < 400/20.9 = 20, resulting in a total icy particle cross section 
of < 0.01 km 2 , which still remains inconsistent with the radar results. 
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Assuming for the moment that the radar cross section properly reflects the total icy 
particle population we compute an upper limit to the total water production rate of < 
2 x 10 27 molec s -1 . The SWAN instrument on the SOHO satellite observes Lya emission 
over fields of view much larger than the radar beam size (10 5 km pixel -1 ). These observations 
are a good point of reference for a total water production rate that is sure to include water 
produced by any large particles observed with Arecibo. |Combi et al. (2011 ) measured a water 



production rate of 6-9 x 10 27 molec s 1 near perihelion (Combi et al.|2011). By this estimate 



it seems that the large particles could account for a substantial fraction of the total water 
production rate of the comet. Note, however, that the SWAN-based water production rates 
are on par with those observed in 4000 km apertures and smaller. Based on the aperture 
sizes and water production rates listed in Table |5j most of the water is produced close to the 
nucleus, perhaps within a few tens or hundreds of kilometers, suggesting that the Arecibo 
observed particles have a low water sublimation rate, if any. 

A dusty large particle population yields a cross section of a = 0.07 — 0.5 km 2 (35J), which 
is in better agreement with the radar results, but suggests that the entire large particle coma 



is within ~ 40 250 km from the comet. There is no indication in Figs. 14 and 17 that the 

large particle coma is truncated on this length scale. Moreover, dusty particles may not have 
sublimating ice. Instead, they could fragment into finer particles. If the large particles are 

dusty, they must fragment on < 40 250 km length scales in order to keep their total cross 

section less than or equal to the observed radar cross section. Just based on the observed 
cross sections, we consider the dusty case to be less likely than the icy case, but still find the 
icy case to be lacking. Of course, a coma of both icy and dusty particles is possible. More 
information on the composition and light scattering properties, including radar wavelengths, 
will be needed to reconcile the radar and MRI observations. 

9. Summary 

Comet Hartley 2 is surrounded by a coma of large particles with radii > 1 cm. With 
observations from Deep Impact, we measured their total flux and flux distribution, based 
on photometry of individual particles. The flux distribution of these particles implies a 
very steep size distribution with power-law slopes ranging from —6.6 to —4.7. We estimate 
that the particles account for 2-14% of the total flux from the near-nucleus coma. The 
spatial distribution of the particles is biased to the anti-sunward direction, as observed by 
the spacecraft both pre- and post-closest approach. Radial expansion from the active areas of 
the rotating nucleus does not explain the observed spatial distribution, even if the ejection 
speeds are very low (~ 1 cm s -1 ). Radiation pressure from sunlight cannot redistribute 
them into the anti-sunward direction on small enough length scales unless the particles have 
extremely low densities (~ 10 -4 g cm -3 ) or low radial ejection velocities (< 10 cm s _1 ). Low 
ejection velocities suggest there should be a strong anti-sunward velocity component in the 
coma, but this does not agree with the velocity distribution observed by |Hermalyn et al. 



(this issue). 



We examined two possible particle compositions. Our models were based on the photo- 
metric properties of the nucleus of Hartley 2 (dusty case: low albedo, 0.3 g cm -3 ) and the 
Jovian satellite Europa (icy case: high albedo, 0.1-1.0 g cm -3 ), and serve as approximate 
limiting cases. 
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The dusty case produces particle size estimates ranging from 10 cm to 2 m in radius, 
the largest of which is just at the limiting resolution of Deep Impact's HRIVIS camera at 
closest approach. Such large particles may be lifted off the nucleus by gas drag if C0 2 is the 
driving gas. Water is a plausible alternative if the water production rate from the nucleus 
is at least 10 27 molec s _1 . Based on the dusty model, the total large particle cross section 
within 20.6 km from the nucleus is 0.07-0.5 km 2 , similar to the 0.89 km 2 radar cross section 



observed by Harmon et al. (2011 ). If these particles are mini-nuclei, we estimate they account 
for 16-80% of the comet's total water production rate (within 20.6 km). However, we can 
all but rule out the dusty case based on total mass estimates of the large particles, which 
are in excess of 10 nucleus masses for densities of 0.3 g cm -3 . Densities < 10~ 3 g cm -3 are 
required to reduce the total mass of the particles to a few percent of the nucleus, which is 
needed to keep the particle mass less than the total mass lost from the comet per orbit (2% 



during the 2010 apparition; Thomas et al. this issue). 



The icy case produces particle size estimates ranging from 1 to 20 cm in radius. These 
particles are easily lifted by water and C0 2 gas drag. Icy particles would sublimate as soon as 
they are heated by sunlight. If the particles have a net sublimation on their sunlit sides, they 
would feel a rocket force that could easily distribute the particles into the anti-sunward direc- 
tion. The sublimation rate excess required for this redistribution is < 10 17 molec cm~ 2 s -1 , 
where the exact value depends on the particle ejection velocities. The cross section of icy 
particles within 20.6 km is much smaller than the observed radar cross section by two to 
three orders of magnitude. The water production rate of the Deep Impact observed particles 
is limited to < 0.1 — 0.5% of the comet's total water production rate. The total mass of 
the particles for a density of 1 g cm -3 is 3-10% the total mass of the nucleus. Thus, porous 
aggregates with p p < 0.1 g cm -3 should be considered more likely than solid ice particles. 

We consider the icy case to be more likely than the dusty case for three reasons: (1) the 
icy particles are more easily lifted by gas drag; (2) we can account for the sunward/anti- 
sunward asymmetry in the particle distribution if ice is sublimating on their sunlit sides; 
and (3) the total large particle mass for the dusty case is much greater than the total mass 
of the nucleus. However, several details are needed in order to test our hypothesis. We need 
an improved icy particle model that treats the large particles more like macroscopic objects 
to better understand their water production rates (if icy) and light scattering properties 
(dusty or icy). We need to better constrain the particle densities, which may rule out the 
dusty case. A hydrodynamic analysis of the near-nucleus coma, especially around the highly 
active small end, would improve our knowledge of the dynamics of large particles. A better 
understanding of the radar coma that includes grains smaller than X/2tc could help resolve 
the discrepancy between our icy particle cross section and the observed radar coma cross 
section. 

If indeed icy with a high albedo, the large particles do not appear to be the source of the 
comet's enhanced water production rate; although, as discussed above, there is much work 
that can be done to refine this conclusion. We suspect that the small icy grains in the jets, 



as observed in Deep Impact IR spectra (A'Hearn et al. 2011, Protopapa et al. 2011), are a 



significant source of water, and the primary cause of the hyperactivity of Comet Hartley 2. 
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Table 1: HRIVIS and MRI images considered in this paper. Table columns are: t — t enci time relative to 
encounter; Exp., exposure time; A, spacecraft-comet distance; Scale, pixel scale at the distance of the comet; 
4>, phase (Sun-comet-spacecraft) angle. 
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Table 2: Online only table. All single power-law flux distributions derived from HR1V1S and MR! 
images. Table columns: F\, the minimum, maximum, and total particle fluxes in the fit; N, number 
of particles in the fit; a, best-fit power-law slope and uncertainty; PkSi the Kolmogorov-Smirnov 
probability that the distribution is a power-law; Ft / F coma , the fraction of the total coma flux 
attributable to point sources with F\ > i*\,min in MRI full-frame data, excluding regions near the 
nucleus. 
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<0.1 
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29.0 
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12.9 
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15.3 
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851 
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425 


-3.08 ±0.11 


41.7 
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4.1 


14.1 
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38.5 
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13.4 
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10.1 
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4.2 


11.4 
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3.9 


16.6 


836 ±7 
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54.5 
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3.8 


18.8 


933 ±7 
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-3.44 ±0.21 


3.0 






mv5004046 


697 


3.6 


12.9 


631 ±5 


124 


-3.50±0.22 


36.4 
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696 


3.0 


13.8 


729 ±6 
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-3.22 ±0.18 


13.5 






mv5004053 


715 


3.1 


21.3 


553 ±5 


110 


-3.36 ±0.22 


99.1 






mv5004056 
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2.8 


11.1 


637 ±5 


135 


-3.34 ±0.20 


35.7 






mv5004058 


807 


3.0 


9.0 


437 ±4 


116 


-3.83 ±0.26 


41.9 
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5.5 


16.4 


1066 ±9 


133 


-4.05 ±0.26 


67.3 


0.0007 


mv5004030 


922 


6.7 


17.8 


884 ±9 


108 


-4.39 ±0.33 


74.5 


0.0006 


mv5004032 


848 


6.0 


22.6 


1773 ±12 


189 


-3.76 ±0.20 


29.8 


0.0012 


mv5004040 


786 


6.1 


24.3 


2053 ±13 


252 


-4.04 ±0.19 


97.8 


0.0019 


mv5004042 


737 


6.0 


24.1 


2483 ± 15 


321 


-4.02 ±0.17 


43.0 


0.0032 


mv5004045 


706 


6.2 


22.2 


2747 ±15 


331 


-3.96 ±0.16 


29.9 


0.0040 


mv6000002 


694 


5.9 


21.1 


3406 ±17 


390 


-3.68 ±0.14 


0.3 


0.0046 


mv5004052 


703 


5.8 


26.2 


3100 ±16 


374 


-3.74 ±0.14 


8.7 


0.0044 


mv5004054 


732 


5.7 


27.1 


2822 ±15 


309 


-3.59 ±0.15 


30.0 


0.0035 


mv5004057 


778 


5.9 


21.4 


1977 ±13 


208 


-3.60 ±0.18 


7.3 


0.0023 


mv6000003 


840 


5.6 


19.8 


2029 ±13 


235 


-3.78 ±0.18 


29.4 


0.0018 


mv5004062 


921 


5.3 


18.3 


1359 ±10 


188 


-4.08 ±0.22 


47.0 


0.0011 


mv5004064 


1004 


4.9 


15.7 


922 ±9 


149 


-3.96 ±0.24 


33.8 


0.0008 










MRI 121 


ins 








mv500402I 


1303 


2.4 


10.6 


2993 ±13 


838 


-3.95 ±0.10 


8.4 


0.0021 


mv5004025 


1107 


2.8 


14.2 


4374 ±16 


1191 


-3.87±0.08 


3.2 


0.0031 


mv5004029 


964 


3.4 


13.9 


6086 ±20 


1291 


-3.88 ±0.08 


0.6 


0.0039 


mv5004031 


885 


3.5 


16.6 


7945 ± 23 


1501 


-3.87±0.07 


<0.1 


0.0048 
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Continued 



Image 
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min 
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mv6000001 


815 


4.1 


19.7 


7337 ±22 


1262 


-3.95 ±0.08 


0.5 


0.0050 


mv5004041 


760 


4.1 


19.8 


9078 ±26 


1468 


-3.91 ±0.08 


0.2 


0.0068 


mv5004044 


719 


3.8 


21.0 


9536 ±26 


1770 


-3.60 ±0.06 


<0.1 


0.0089 


mv5004046 


697 


4.3 


25.4 


9012 ±25 


1391 


-3.76 ±0.07 


0.1 


0.0087 


mv5004051 


696 


3.7 


29.2 


9695 ±26 


1769 


-3.65 ±0.06 


<0.1 


0.0104 


mv5004053 


715 


3.5 


23.1 


9603 ±26 


1769 


-3.57 ±0.06 


<0.1 


0.0101 


mv5004056 


753 


3.4 


20.7 


9543 ±26 


1806 


-3.58 ±0.06 


<0.1 


0.0083 


mv5004058 


807 


3.2 


16.8 


10378 ±27 


1937 


-3.65 ±0.06 


<0.1 


0.0073 


mv5004061 


882 


3.0 


24.0 


8495 ± 24 


1681 


-3.65 ±0.06 


<0.1 


0.0056 


mv5004063 


962 


2.6 


14.3 


6967 ±21 


1771 


-3.64 ±0.06 


1.2 


0.0047 


mv5004066 


1096 


2.4 


14.7 


4796 ±17 


1426 


-3.82 ±0.07 


1.2 


0.0034 


mv5006001 


1294 


2.1 


12.2 


2373 ±11 


965 


-3.87±0.09 


17.8 


0.0021 










MRI 501 


ms 








mv5002051 


5148 


0.4 


1.1 


54 ±1 


105 


-4.54 ±0.35 


82.3 


0.0002 


mv5002061 


4082 


0.5 


1.7 


102 ±1 


192 


-4.70 ±0.27 


92.6 


0.0003 


mv5002069 


3210 


0.7 


3.0 


243 ±2 


292 


-4.26 ±0.19 


5.8 


0.0005 


mv5004001 


2523 


0.8 


3.4 


653 ±3 


721 


-3.76 ±0.10 


11.6 


0.0011 


mv5004005 


2172 


0.9 


5.6 


1052 ±5 


837 


-3.70 ±0.09 


17.6 


0.0013 


mv5004009 


1858 


1.1 


6.7 


1526 ±7 


1034 


-3.69 ±0.08 


0.1 


0.0017 


mv5006016 


1886 


0.9 


7.3 


1288 ±6 


991 


-3.59 ±0.08 


9.0 


0.0015 


mv5006020 


2219 


0.7 


4.3 


1028 ±5 


954 


-3.72 ±0.09 


4.2 


0.0013 


mv5006024 


2643 


0.7 


3.2 


400 ±3 


411 


-3.81 ±0.14 


39.6 


0.0007 


mv5006030 


3345 


0.7 


2.8 


166 ±2 


184 


-4.04 ±0.22 


63.0 


0.0003 


mv5006037 


4187 


0.5 


1.6 


81 ±1 


120 


-4.85 ±0.35 


37.5 


0.0002 



a The absolute calibration uncertainties are 5% for HRIVIS, and 10% for MRI. 
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Table 3: Summary of single power-law fits to individual HRIVIS and MR! flux distributions listed in Table [2] 
Table columns: F\, mean minimum and maximum particle fluxes; N, mean number of particles; a, minimum, 
maximum, mean, and standard deviation of the flux distribution power-law slopes; PkSi mean Kolmogorov- 
Smirnov probability that the distributions are power-laws; Ft / F comai mean fraction of the total coma flux 
attributable to point sources with F\ > F\_ min in MRI full-frame data, excluding regions near the nucleus. 



Instrument 



A 



Fx 



N 



mm max mm max 
(km) (icr 12 



max mean 



Pks FT/F coma 
(%) 



WirT 



U1T1 



- 1 ) 



HRI 

MRI sub-frame 
MRI 41 ms 
MRI 121 ms 
MRI 501 ms 



694 
696 
694 
696 



915 

964 

1006 

1303 



1.3 
3.6 

5.5 

3.1 



30.0 
14.2 
20.2 
17.0 



545 

127 

245 

1489 



-3.08 -2.67 

-4.31 -3.22 

-4.39 -3.59 

-3.95 -3.57 



1858 5148 0.7 



3.5 531 -4.85 -3.59 



MRI (distant) a 2172 2643 9.0 46.9 731 -3.81 -3.70 



-2.87 
-3.70 
-3.90 
-3.76 
-4.06 
-3.75 



0.14 
0.35 
0.22 
0.14 
0.43 
0.04 



0.125 
0.341 
0.384 
0.021 
0.330 
0.182 



0.0023 
0.0059 
0.0008 
0.0011 



a Averages of the the four 501 ms MRI images with A > 2000 km and TV > 300 particles, with which we 
use to define the fraction of the coma flux attributable to large particles. For this row only, F\ has been 
corrected to a distance of A = 700 km before averaging. 



Table 4: Flux range (-F^mim ^A max); radius range (ai cy ), total number (N), flux (-Fa,t); fraction of coma 
flux in particles (Ft / 'F coma ) , cross section {(Ti C y), mass (Mi cy ), and water production rate (Q(H20)) for the 
large particles measured at A = 2100 — 2600. The radius, cross section, mass, and water production rate 
are computed using an icy composition with Europa's photometric parameters (at <p = 86°) and a density 
of 1.0 g cm~ 3 . We also extrapolate these quantities down to 10 -14 W m -2 urn -1 , and to our total coma 
measurement for a 20.6 km aperture at A = 5148 km. Uncertainties are derived from the power-law slope 
uncertainties (±0.02) and are symmetric in log space. 



r? a j? a 

^A.min r A,max 

xicr 12 

(W m~ 2 jim -1 ) 



&icy 

(cm) 



loglO N lo §10 ^A,T a 

(W m~ 2 urn -1 ) 



(%) 



loglO ff icy 
(cm 2 ) 



l°glO M icy 
(g) 



1°§10 Wmax 

(s-i) 



A m 2400 km, full-frame (excluding masked regions) 



9.00 b 


45 b 


7.4 


1.00 


45 


2.5 


0.10 ° 


15 


0.8 


0.01 


45 


0.2 



16.5 2.86 b -8.06 b ' c 

16.5 4.90±0.02 -6.86±0.02 

16.5 6.76±0.04 -5.97±0.03 

16.5 8.61±0.06 -5.11±0.05 



0.11 b 

1.8±0.1 
13.8±1.1 

98.7 



+ 12.8 
11.3 



5.22 9.07 22.56 

6.42±0.02 9.92±0.01 23.77±0.02 

7.32±0.03 10.39±0.03 24.66±0.03 

8.17±0.05 10.79±0.04 25.52±0.05 



A = 5148 km, 20.6 km aperture 



1.00 
0.10 ■ 



15 
15 



2.5 
0.8 



16.5 
16.5 



6.28±0.03 
9.02±0.07 



-6.69±0.02 
-5.79±0.03 



1.8±0.1 

13.8±1.1 



6.59±0.02 
7.49±0.03 



11.30±0.03 

12.66±0.06 



23.94±0.02 
24.83±0.03 



a Fluxes have been corrected to a distance of 700 km. 

b Measured values. 

c The instrumental and calibration uncertainties are 2 X 10 -12 W m -2 urn -1 and 10%, respectively. 

d Assumed value. 

e Our preferred lower flux limit is between 0.1 and l.OxlO - 12 W m -2 urn -1 . 
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Table 5: Summary of comet Hartley 2 water production rates observed near perihelion, and derived in this 
study. 



Date(s) Aper. a 

(UT, 2010) (km) 



Q(H 2 0) 
'10 27 molec s _1 " 



Notes 



18 Oct 

19 Oct 
22 Oct 
23-31 Oct 
31 Oct b 

1 Nov 

2 Nov 
4 Nov 
4 Nov 



9.1 x 10 5 

19 

19 

1.6 x 10 4 

3800 

1.1 x 10 6 

300 

22 

20.6 



8.70 ±0.38 
3.45 ±0.04 
6.78 ±0.26 
<2. 

11.6±0.7 
6.38 ±0.12 
10. 
11.45 ±0.65 
<0.004 



Combi et al. 



2011 



Mumma et al. 



Mumma et al. 



2011 



2011 



Large particles only, 



Harmon et al. 



Knight and Schleicher this issue 
Combi et al."||201l ' 



2011 



+ this work 



A'Hearn et al.||2011 



Dello Russo et al. 



Large particles only, Table |4| 



2011 



Slit half-width or aperture radius. We have assumed the 0.43" slit for the Mumma et al. (2011 1 observations, 



and a 3 pixel radius for the Combi et al. (2011) observations. 

b Mean and standard deviation of 5 measurements taken on 31 Oct. 
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